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ABSTRACT 


It  has  been  suggested  by  the  author  that  some  aspects  of 
severely  distorted  flow  into  multistage  compressors  may  be  exam- 
ined utilizing  an  integral  technique.  The  general  idea  of  the  proposed 
technique  is  clear  enough;  the  appropriate  equations  of  motion  and 
energy  are  integrated  peripherally  and  radially,  using  reasonable 
assumptions  for  the  distributions  of  velocity  and  thermodynamic 
properties,  and  thereby  reduced  to  ordinary  non-linear  differential 
equations  for  the  parameters  that  describe  the  distributions.  The 
questions  that  arise  are  whether  the  cascade  characteristics  may 
be  described  appropriately  over  wide  variations  of  inlet  angle,  in- 
cluding stall,  and  whether  the  profiles  may  be  characterized  by  a 
sufficiently  small  number  of  parameters  to  make  the  technique 
attractive. 

The  present  paper  examines  a specific  example  of  distorted 
inlet  flow  through  the  two-dimensional  annulus  of  a multistage  com- 
pressor which  can  be  solved  completely.  It  is  shown  that  the  essen- 
tial features  of  this  exact  solution,  including  stall,  may  be  described 
by  a two-parameter  family  of  profiles  and  that  an  integral  technique, 
utilizing  these  elementary  profiles,  will  yield  essentially  the  same 
results.  While  it  is  not  clear  that  comparable  success  would  hold  for 
the  three-dimensional  problem,  the  results  confirm  the  contention 
that  the  two-dimensional  problem  may  be  treated  with  acceptable 
accuracy  by  an  integral  technique. 


1.  INTRODUCTION 

Non-uniform  flow  conditions  entering  multistage  compres- 
sors are  among  the  most  significant  factors  that  limit  the  perform- 
ance of  these  units  and  restrict  the  operation  of  engines  of  which 
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they  are  components.  ’ Under  some  circumstances  the  non- 
uniformities are  attenuated  in  the  first  few  stages  while  otherwise 
they  may  persist  through  the  entire  machine,  causing  severe  per- 
formance degradation  and  unacceptable  blade  loads.  The  develop- 
ment or  smoothing  out  of  severely  distorted  inlet  states  is  a complex 
cooperative  phenomenon  among  many  successive  stages  ’ ’ and  is 

not  readily  understood  in  terms  of  individual  blade  row  performance. 
It  is  the  aim  of  analytical  work  in  this  field  to  formulate  such  pro- 
lems  within  the  requisite  accuracy  and  to  find  methods  by  which 
physically  relevant  results  may  be  extracted.  It  is  the  contention 
of  some  workers,  including  the  author,  that,  preferably,  such 
formulations  and  methods  should  lead  relatively  simply  from  a 
physical  model  to  significant  results. 

One  contender  for  the  appropriate  analytical  technique  is 
an  integral  method  in  which  the  equations  of  motion  and  energy  are 
integrated  over  the  compressor  annulus  at  a fixed  point  along  the 
compressor  axis.  The  most  familiar  example  of  this  technique  is 
the  Karman  integral  method  for  the  boundary  layer,  although  the 
general  idea  has  found  wide  application.  The  applicability  of  the 
technique  to  the  problem  of  asymmetric  compressor  flow  rests 
not  only  upon  whether  it  is  simpler  than  a detailed  numerical  cal- 
culation, but  even  more  critically,  whether  the  velocity,  pressures 
and  temperature  profiles  may  be  described  adequately  by  a small 
number  of  parameters  and  whether  the  blade  characteristic  can  be 
specified  reasonably  over  a wide  range  in  angle  of  attack.  One  can 
not  settle  these  questions  for  all  circumstances,  but  one  example 
in  which  an  "exact"  calculation  and  an  integral  method  are  in 
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agreement  would  be  a guide  to  the  merits  and  limitations  of  the 
method. 

It  is  the  aim  of  this  paper  to  construct  such  an  example  in 
which  the  exact  solution  can  be  carried  out.  The  problem  utilizes 
a blade  row  and  stalling  characteristic  and  a solution  technique 
introduced  by  the  author^  ^ several  years  ago  in  the  analysis  of 
propagating  stall  in  single  blade  rows.  To  make  matters  simpler, 
the  example  considers  a two-dimensional  annular  section  of  the  com 
pressor,  represented  as  a repetitive  flow  pattern  in  the  plane,  and 
the  fluid  field  is  regarded  as  uncompressible  and  inviscid  although 
the  individual  blade  rows  impose  losses  on  the  flow  across  them. 

The  essential  questions  to  be  asked  of  the  example  is  whether  the 
profiles-especially  the  velocity  profile-can  be  described  in  terms 
of  a small  number  of  parameters.  If  the  required  number  is  few, 
then  the  possibility  exists  that  the  integral  approach  may  prove 
useful  in  treating  problems  of  distorted  inlet  flow. 
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2.  BLADE  CHARACTERISTICS  AND 
REPRESENTATION  OF  STALL 

Consider  the  flow  of  an  incompressible  fluid  through  a 
two-dimensional  cascade  utilizing  the  notation  of  Figure  1;  (3,  and 
(3Aare  the  relative  inflow  and  discharge  angles  under  conditions  of 
steady,  uniform  flow  and  U is  the  constant  axial  velocity.  We  re- 
place the  cascade  by  an  actuator  surface  coincident  with  the  axis. 
Under  these  circumstances  the  pressures  ahead  and  behind  the  cas- 
cade are  uniform  and  denoted  -f>,  and  7°i.  respectively.  The  charac- 
teristics of  the  cascade  will  be  described  by  the  turning  angle, 

(3,.-  (1,  , and  the  static  pressure  rise,  TV  "P.  £ &V , in  terms  of  the 
inlet  angle  (3,  . We  shall  take  these  characteristics  to  be  those 
given  in  Figures  2 and  3 , a form  introduced  in  Ref.  6 , for  the 
analysis  of  propagating  stall.  The  pressure  rise,  , is  a discon- 

tinuous function  of  inlet  angle  and  drops  to  zero  at  the  local  inlet 
stalling  angle,  |3*,  but  is  continuous  and  smooth  for  angles  <■  (3* 
The  outlet  angle,  , is  a smooth  function  of  p,  through  the  stall. 
The  stall  is  thus  characterized  by  the  absence  of  static  pressure  rise 
in  the  separated  blade  channel,  that  is,  the  relative  diffusion  drops 
to  zero. 

Now  we  consider  non-steady  flows  in  the  cascade  to  be  peri- 
odic over  a range  xit  R , corresponding  to  the  periodicity  about  the 
compressor  annulus,  and  to  be  reducible  to  steady  flow  by  an  ap- 
propriate velocity  of  translation  along  the  y-axis.  A disturbance 
that  enters  from  a fixed  obstruction  ahead  of  the  inlet  is  steady  with 
respect  to  the  stationary  blade  rows,  but  we  must  move  with  a verti- 
cal velocity  ft  with  respect  to  a rotor  blade  row  to  achieve  a 
steady  flow  field.  In  the  case  of  self-induced  disturbances,  such  as 
a stall  propogation  that  moves  with  an  angular  velocity  u>  , we  must 
move  with  a velocity  - con  along  a stator  and  with  a velocity 
along  a rotor  in  order  to  achieve  a steady  flow  field.  This  is 
possible  because  the  actuator  surface  concept  suppresses  effects 
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oi  individual  blades  and  these  are  the  only  identifiable  feature's  of 
blade  row  geometry. 

Where  the  appropriate  transformation  along  the  ^ -axis 
has  been  made,  we  shall  call  the  angles  of  the  steady  flow,  <9,  and 
&z.  upstream  and  downstream  respectively;  the  relations  between  the 
relative  flow  angles  ^ and  the  steady  absolute  angles  0,  t 

0^  are 

j-  _ , t*JR  1 

Cwvi  &,  - {},  - y 

for  a stator  and 

/-  a r a - L 

LanG,-  Canfi,  - - -j- 

for  a rotor  where  oo  is  the  absolute  angular  velocity  of  the  disturbance 
and  -O.  is  the  angular  velocity  of  a rotor.  For  a stationary  disturbance 
introduced  at  the  inlet,  we  set  wo  ; for  a propagating  stall,  we  de- 
termine oo  as  a characteristic  value. 


A 
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3.  NON-UNIFORM  FLOW  IN 
A SINGLE  BLADE  ROW 

Under  any  of  the  circumstances  just  considered,  the  flow 
field  induced  by  a single  blade  row  may  be  described  as  in  Figure  4, 
in  which  any  existing  stalled  region  of  length  (xWG)*-  is  located 
symmetrically  with  respect  to  the  * -axis.  The  flow  angles  <9,  and 
^ and  the  tangential  velocities  and  are  those  that  exist 

in  the  absence  of  stall.  We  are  concerned  with  small  deviations  from 
these  two  uniform  fields  introduced  by  a stall  or  an  upstream  dis- 
tortion. 


The  perturbation  equations  applicable  to  either  upstream 
or  downstream  flow  fields,  are 


du.  ^ O'T 
ax  d + 
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ax  v 
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The  disturbance  S'  of  the  flow  angle  from  the  steady  flow  angle 
may  be  expressed  in  terms  of  the  velocity  perturbations  and  if  we 
retain  only  terms  linear  in  perturbation  quantities,  this  relation  is 


(i-f  ' 


y 

U 


du-  . 

Now  if  we  utilize  equation  3 to  eliminate  the  term  &&  in  favor  ot 
in  the  x-component  of  momentum  conservation,  equation  4, 
we  obtain  - (/^  (-jj  - £ j ° and,  upon  introduc- 

ing the  flow  angle  perturbation,  equation  6, 


0* 


Cmt*  & ) t/ 
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Similar  operations  on  the  second  equation  of  motion,  equation  5, 
yields 

, , . id#  ^ 

(t  * ° a x ' ? 
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If,  therefore,  we  choose  two  new  dependent  variables 

& ' J 


P 


10 


the  two  equations  of  motion,  in  the  forms  of  equations  7 and  8, 
become 


aP 

- 
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5 X 

aP 

d® 
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c>  > 

a * 

Thus  the  complex  variable  p*>  & 

may  be  considered  an  analytic 

function  of  the  comple: 

x variable  x-> 

"i  ■ * 

0 ' Wcf) 
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This  description  of  the  pressure  and  flow  angle  perturbations 
in  terms  of  an  analytic  function  does  not,  of  course,  necessitate  that 
the  flow  be  irrotational.  The  irrotational  part  of  the  field  is  deter- 
mined by  U/<2',  but  additional  "shear  flow"  solutions  may  be  super- 
imposed because,  when  the  flow  is  steady,  such  flows  carry  no  pres- 
sure field  and  induce  no  angle  variations.  In  the  present  problem 
the  admissible  shear  flows  are  easily  obtained  by  writing  a linear 
combination  of  equations  4 and  5 to  give 


cos©  — 
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where  the  operator  denotes  differentiation  in  the  streamwise  direc- 
tion. The  quantity  in  brackets,  which  is  proportional  to  the  pertur- 
bation of  stagnation  pressure,  therefore  varies  only  from  one  stream- 
line to  the  other,  and  the  integral  of  equation  14  may  be  written 


t-<w>©  -*  - 
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From  equations  15  and  6 it  is  evident  that  t is  the  vorticity  and 
that  it  is  transported  along  the  undisturbed  streamlines.  In  fact, 
from  these  same  equations,  the  two  velocity  components  may  be 
written 


C*  ~ & & 


(I  -* 
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Now  let  us  consider  a single  blade  row,  represented  by  an 
actuator  surface  located  on  the  ^-axis,  for  which  the  undisturbed 
uniform  flow  moves  at  an  angle  0,  upstream  and  0x.  downstream 
in  accordance  with  the  diagram  of  Figure  4.  There  is  a corres- 
ponding pressure  rise  A) o across  the  actuator.  Let  us  suppose 
that,  due  eithe  r to  a self  induced  or  upstream  induced  disturbance, 
the  cascade  stalls  over  a fraction  oc  of  its  circumference,  and  that 
the  stall  region  is  located  as  shown  in  Figure  4.  According  to  the 
chosen  stall  model,  the  stall  induces  no  disturbance  to  the  relative 
discharge  angle  from  the  cascade,  but  the  pressure  rise  Ap  drops 
to  (say)  zero  across  the  stalled  portion.  The  condition  on  the  pres- 
sure field  is  then  that 

£ I-*  ) ($“*  (0,  g ) - f L*  ( 1+tflL"1©,  1 t »,  j ) 

[ - Ay*  ; n-  4 ‘ -2-  *«■»■#• 

1 *■»■««.  j y 

so  that  the  difference  between  the  solutions  for  the  pressure  fields 
upstream  and  downstream  consists  of  positive  jumps  at  the  points 
^ . (*,♦•£)  kv R and  equal  negative  jumps  at  the  points  (*-£)**& 

Since  6^  and  are  harmonic  functions,  one  may  begin  to  con- 

struct a solution  by  obtaining  the  harmonic  function  that  represents 
the  value  of  Q <•».  jp  ' an  1 is  regular  elsewhere.  Now 
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the  imaginary  part  of  the  function  £ log  has  a positive 

unit  jump  at  for  X-  °*  and  a corresponding  negative  jump 

for  X - o-  . The  sort  of  function  we  require  may  then  be  con- 
structed by  situating  the  appropriate  logarithmic  singularities  at 

and  the  negative  of  such  functions  at  If- * tn-Siio® 
The  analytic  function 


1 i 1*"  n [ jut  n ' *•  J I 

c r 12)  • v Joj  — - ~~~  ( 

* l n 1 ■*  *’  * 3 ) 

posesses  this  property,  since  the  hyperbolic  sine  in  the  numerator 
has  simple  zeros  on  the  j -axis  at  x • Figure  5,  and  the 

denominator  has  corresponding  simple  zeros  on  the  g -axis  at 
JL  «■  vw  — . The  values  of  the  real  and  imaginary  parts 

JLtWX  *■ 

along  the  ^ -axis,  which  are  involved  in  satisfying  matching  condi- 
tions across  the  actuator  line,  are  shown  also  in  Figure  5.  The 
real  part,  F<x,  ft  is  single  valued  but  exhibits  the  expected 
logarithmic  singularities  at  n * • The  imaginary 

part,  Go,)!  with  the  branch  cut  taken  along  the  "j  -axis,  exhibits 
a jump  of  magnitude  2 in  the  intervals  n-  a £ -4-^  t a and  is  con- 
tinuous elsewhere.  Far  upstream  of  the  cascade,  as  X— * - $ 

6c-oo  has  a constant  value  oc  and  far  downstream,  as 

has  a constant  value  - <x  . 

Now  the  solutions  (?,u.p  a.r\JL  (?*.<.*.  j\  for  the  pres- 
sure field  can  be  constructed  from  a multiple  of  G < *,  J ' and  an" 
other  harmonic  function  which  must  be  regular  except  at  the  singu- 
larities along  the  imaginary  axis.  Hence  the  supplementary  func- 
tion can  differ  from  a multiple  of  J ) by  at  most  a constant. 

The  flow  angle  and  pressure  far  ahead  of  the  cascade  can  vary  by 
only  constant  magnitude  from  the  undisturbed  ; these  values  corres- 
pond to  variations  in  flow  angle  and  pressure  imposed  far  upstream. 
If  we  assume  these  to  vanish,  then  the  prsssure  and  angle  pertur- 
bation fields  become 
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where  the  constants  A , B , Z>  remain  to  be  determined.  The  avail- 
able matching  conditions  are  (i)  that  the  axial  velocity  is  continuous 
across  the  cascade  and  (ii)  that  the  relative  discharge  angle  is  given 
in  terms  of  the  inlet  angle  according  to  Figure  3.  For  the  present 
consideration  we  shall  assume  that  the  relative  discharge  angle  is 
fixed,  independent  of  the  inlet  angle.  The  representations  of  the 
axial  velocities  upstream  and  downstream  of  the  cascade  may  be 
written 

r ~ Pt  ~ ®,  to^e,  23 

XT  - - («l  * - (*- 

■*  a-t.  F < *>.>-*  + lo*.  [ 6 (o«,  7.--*  t««a  ) ♦ * ] 24 


Where  u-thas  been  supplied  with  an  additive  constant  which  assures 
that  u averages  to  zero  far  downstream.  The  condition  that  the 
relative  discharge  angle  from  the  cascade  be  invarient  requires  ex- 
pressing the  relative  flow  angle  perturbation  in  terms  of  perturba- 
tions seen  from  the  steady  flow  field.  Referring  to  equations  1 and  2, 
which  relate  these  flow  angles. 
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so  that  to  the  first  order  in  small  disturbances, 

r /wp  i \jJ 


/wf  \ u- 

l U / 
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The  condition  that  the  outlet  angle  be  unchanged  is  therefore 
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which  we  write  in  the  form 
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Now  all  of  the  field  quantities  are  represented  by  equations  19 
through  24  and  substitution  into  the  two  matching  conditions  provi- 
des relations  among  the  unknown  coefficients.  In  fact,  the  coeffic- 
ients of  Fio.p  and  [Gc*.*)-**  ] must  vanish  identically  because 
these  two  functions  are  orthogonal  over  the  range  o s - -if'? 
Substituting  into  the  expression  -jj1-  < \)  I r {7*  to'3'  for  continuity 

of  axial  velocity  leads  to  the  algebraic  equations 

I ~ , 1 \ / hn  G,  ^ \ — _ _ _ 3 0 


r 


lirnie.  7^‘»J-4Uu^  i+t'-'sj6’  b-  i+iS* 

Similarly,  substitution  into  the  condition  on  the  cascade  discharge, 
angle,  equation  28,  gives  the  second  pair  of  algebraic  equations. 

_L*_  A «.  * Kai  - 32 

utoJ'e.  h /•*  tm~  ?u 
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These  equations  constitute  a linear  set  for  A,  B,  «*». , fa,.  and  may 

be  written,  after  some  manipulation,  in  the  matrix  form 

fi  i i ^ \ If  1 f +•»•»  1 
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The  determinant  of  the  matrix  on  the  left. 


f « \\  . 

/ kus,  • V 

l 1 ■*  tavfa,  \ ■*  t*.1©.  / 

is  non-vanishing  and  the  resulting  values  of  AJ  B,  <*■» , b^ 
proportional  to  , the  pressure  loss  coefficient  at  stall. 

The  results  will  not  be  written  down  explicitly  now  because  they 
will  emerge  later  in  a more  general  form.  It  is  well  to  note. 


-13- 


k— 

i 


however,  that  because  a.u*o  , loj.  / o , the  flow  far  down- 
stream of  the  cascade  is  permanently  distorted  and  corresponds  to 
a wake  generated  by  the  vorticity  sheet  from  the  stalled  zone. 

The  physical  problem  of  the  single  blade  row  is  still  in- 
complete. It  has  been  assumed  that  the  cascade  is  stalled  over  a 
certain  portion  of  its  periphery  but  the  conditions  which  lead  to  the 
stall  have  not  been  considered.  We  must  examine  two  possibilities: 
i)  the  stall  is  induced  by  a flow  distortion  originating  upstream  of 
the  cascade  and  ii)  the  stall  is  self-induced.  The  case  of  the  self- 
induced  or  "propagating  stall"  will  be  presented  in  this  section;  the 
stall  induced  by  flow  distortions  will  be  treated  subsequently. 

Now  whether  the  cascade  stalls  locally  is  determined, 
according  to  Figure  2,  by  the  local  inlet  angle  (3,  . In  the  absence 

of  disturbance  the  uniform  inlet  angle  is  assumed  less  than  the 

stall  angle  and  in  the  presence  of  the  stall,  the  perturbation 

^3  must  relate  appropriately  to  the  angle  increment  (3,*- 

required  to  stall  the  airfoil.  Using  the  relationship  we  have  found 
between  absolute  and  relative  perturbation  angles,  we  have 


0*  i^,  g),  f*,j)  + + h£  XT’1®'*' 

Now  using  the  solution  forms  we  have  found  for  <9 , (x. gj  • ( x.  u l 
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The  stable,  self-induced  stall  cell  requires  that  the  inlet  flow  angle 
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change  from  a value  below  p*  to  a value  equal  to  or  greater  than 
(3,  * at  the  edges  of  the  cell  and  this  condition  can  be  satisfied  only 
if  the  coefficient  of  vanishes.  Thus 

( id»«l6i)KA  ■*  -(ntJ©4xiaue,j  e ;o 


37 


This  additional  condition  on  the  constants  A and  & may  be  taken 
together  with  the  fourth  of  equation  34,  another  homogeneous  condi- 
tion, to  require  that,  for  a non-trivial  solution. 
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which  constitutes  a characteristic  value  condition  for  K . Physically, 
this  states  that  the  self-induced  stall  cell  must  move  with  a certain 
speed  in  order  that  the  flow  field  satisfy  the  condition  on  the  inlet 
angle  approaching  the  stall  cell.  This  free  or  self-induced  stall 
cell  has  become  known  as  propagating  stall  because  it  must  move 
along  the  blade  row  at  this  speed  to  assure  the  stability  of  the  lead- 
ing trailing  edges  of  the  stall  cell.  The  value  of  the  determinant  in 
equation  38  gives  the  characteristic  value  as 
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where  the  relation  between  relative  and  absolute  flow  angles  has  been 
used  to  obtain  the  final  expression.  Here  ^3,  is  essentially  the  rela- 
tive inlet  stalling  angle  of  the  blade  row  in  question  and  coH  is  the 
velocity  of  stall  propagation. 


Note  that  the  value  of  K determined  for  the  self-induced 
stall  does  not  make  the  matrix  in  equation  34  singular  but  allows 
determination  of  the  complete  set  of  coefficients,  in  particular, 
the  values  of  and  l»xwhich  give  the  stall  wake  that  appears  far 
downstream  from  the  blade  row. 
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4.  SINGLE  BLADE  ROW;  GENERAL 
ANALYSIS 

The  analysis  of  section  3 has  shown  that  the  vorticity  shed 
from  a stalled  blade  row  produces  a wake  that  deforms  the  flow  field 
far  downstream  of  that  blade  row.  It  follows,  moreover,  that  sub- 
sequent blade  rows,  in  relative  motion  to  the  stalled  region,  will 
experience  strong  deviations  of  the  relative  inlet  angle  which  may 
induce  stall  in  that  blade  row  also.  In  general,  then,  solution  of  the 
multistage  problem  will  require  knowledge  of  the  behavior  of  single 
stages  with  adeformed  input.  It  is  the  aim  of  this  section  to  find  such 
a solution,  suitable  for  use  in  multistage  analysis. 

In  the  previous  section  it  appeared  that  the  axial  velocity 
profile  from  a stalled  single  blade  row,  consisted  of  a linear  com- 
bination of  Fto,  and  &(o*.  y. -r  +*  even  though  the  flow 

angle  perturbation  vanished  far  downstream.  Therefore  a distor- 
tion of  this  sort  will  appear  upstream  of  the  next  blade  row.  In  gen- 
eral, then,  a blade  row  will  have  some  linear  combination  of 

and  6(o-t  frta-9)-*  upstream;  a solution  that  accommo- 
dates this  initial  condition  will  permit  also  a distortion  of  this  type 
to  be  prescribed  ahead  of  the  entire  compressor. 

To  satisfy  these  requirements,  the  pressure  and  angle  of 
perturbations  may  be  written 
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while  the  corresponding  axial  velocity  distributions  are 

77  E ‘ " 77 TL'e,*}  ~ ~77Z-'o?)  7~*6’- 

+ o>,  F(o,}-*tc~if.)  - t,  [0(0-,^- 1 -<*  ) 
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!±i  : - /p  „ AX  - c \ _ / A „ \ 

U vu*  fu*  HttvJSL  /^«.|\  j > 4 

■*  o.j  F ( o,  1J  -*  fro.*  + b^G  C®«,  •j-xtw-  ©,T  +<x  7 

The  prescribed  conditions  ahead  of  the  blade  row  determine  the 
values  of  a,,  fe, } C , J>,  and  oc  and  the  values  of  A , B , a* 

andA  are  required.  These  constants  are  evaluated  as  a 
consequence  of  matching  the  axial  velocity  perturbations. 


to.  -jS 

and  fixing  the  relative  discharge  angle 

@x  lo,'j\  ■*  K =o  47 

The  pressure  loss  caused  by  stall  has  already  been  accounted  for 
in  constructing  the  solutions.  Now  a prescribed  inlet  distortion, 
fixed  by  the  values  of  *.  and  lo,  , may  or  may  not  induce  stall 
over  a portion  xnnec  of  the  blade  row  periphery.  This  ques- 

tion is  determined  by  the  relative  inlet  angle;  if  stall  is  indicated, 

Ay,  in  equations  40-45  is  set  equal  to  the  value  appropriate  for 
the  blade  row,  if  not,  a-jo  is  set  equal  to  zero. 

Another  feature  that  appears  when  the  stall  is  forced  by  the 
inlet  distortion  is  that  the  value  of  K is  prescribed.  For  a station- 
ary distortion,  such  as  strut  or  a separated  inlet,  c • o , but  in 
the  unusual  circumstances  of  a moving  disturbance,  which  may  occur 
with  an  ingested  vortex,  the  value  of  aj  is  assumed  known.  This 
will  mean  that,  in  general,  the  inlet  angle  distortion. 
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will  not  be  a simple  square  wave  but  will  contain  the  singu- 
larities characterizing  j)  • The  condition  for  stall  in  this 

circumstance  will  be  discussed  later. 

Substitution  into  conditions  46  and  47  leads,  in  a straight- 
forward manner,  to  the  result  that  -Z^_  •=•  o and  a set  of  linear 
algebraic  equations  for  the  unknown  quantities  A,  B,  *•«.,  bk. 

These  are  identical  with  the  equations  34  except  for  additional  terms 
on  the  right-hand  side 
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where  and y u are  functions  of  O, , , K and  are  tabulated 

in  the  Appendix. 

Now  a,  and  b,  are  given  through  prescription  of  the 
flow  far  upstream  but  yjz  is  either  a known  value  or  zero  accord- 
ing to  the  approach  angle  relative  to  the  blade  row.  The  relative 
inlet  angle  perturbation,  given  by  equation  36,  becomes  for  the 
present  example. 
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the  average,  contribute  to  the  inlet  angle  in  the  region  under  con- 
sideration. The  function  G (°">  J 1 - * however,  equals  /-*  in  the 
stalled  region  and  - « in  the  region  exterior  to  the  stall.  In  detail, 
then,  the  inlet  angle  perturbation  is,  in  the  distorted  region 
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where  A is  the  determinant  given  in  the  Appendix.  Then  if 

- p.'-/3-  52 

the  cascade  will  stall  in  the  region  ~tja«  ~ - n n‘K  in  which  the 

upstream  distortion  impinges  upon  it  and  the  solution  is  correct  as 


it  stands.  If,  however,  the  inequality  52  is  not  satisfied  then 
the  solution  given  by  equations  49-52  must  be  evaluated  with 

set  to  zero.  This  condition  also  determines  the  pressure 
loss  far  downstream  from  the  cascade,  according  to  equation  42, 
as  -^T*  x or  zero. 


5.  ANALYSIS  OF  MULTISTAGE  COMPRESSOR 
WITH  WIDE  SPREAD  BLADE  ROWS 


The  results  of  the  last  two  sections  demonstrate  that  a 
compressor  blade  row,  represented  by  an  actuator  surface,  which 
is  subjected  to  a distortion  in  axial  velocity  profile  that  is  a linear 
combination  of  and  G(o-,p-ot  , produces  a distorted  dis- 

charge flow  which  is  another  linear  combination  of  the  same  two 
functions.  This  general  result  holds,  whether  or  not  the  cascade 
in  question  stalls  as  a result  of  the  distorted  input.  It  follows, 
then,  that  we  can  construct  the  flow  in  a multistage  compressor, 
represented  in  this  two-dimensional  approximation,  from  linear 
combinations  of  the  two  fundamental  profiles  given  above.  The 
situation  holds  regardless  of  how  many  of  the  various  blade  rows 
stall. 


Number  the  blade  rows  from  the  inlet,  the  first  of  which 

may  be  the  inlet  vanes,  the  successive  ones  may  alternate  rotor 

th 

and  stator  rows.  Number  the  region  ahead  of  the  n blade  row 
by  the  subscript  n , that  is  «•■»,  (jJ,  P„_6l,and  the  region  downstream 
by  the  subscript  n*i  , that  is  “■•Ml  > . 

Associated  with  each  of  these  blade  rows  is  a stalling  inlet  angle 
|3„*  , a fixed  relative  outlet  angle  (3n+,  > a pressure  loss  at 

stall,  A-jo„  , and  a value  K„  of  the  blade  motion  parameter.  For 
flow  fields  generated  by  a stationary  upstream  disturbance,  the 
values  of  are 


- A n/o 
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Rotor. 
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1 -t  t CMA*-  0J, 

We  presume  the  unperturbed  state  to  be  known  and  this  means,  in 
detail,  that  values  of  U,  , (3„  are  known.  The  perturbations 

we  wish  to  compute  are  caused  by  a distortion  upstream  of  the  com- 
pressor inlet  described  by  the  values  of  CL,  and  b,  ; generally 


this  takes  the  form  of  a distortion  to  the  axial  velocity  profile  given 
by 


b,  [ 6 c 3 j - ] 


We  shall  assume,  for  the  purpose  of  this  demonstration,  that  the 
pressure  and  angle  disturbances  vanish  far  upstream  and  that 

. Note  specifically  that  there  is  no  perturbation  to  the  total 
flow  rate,  that  is 


u,  otj  c 


so  that  the  mean  or  reference  axial  velocity  is  undisturbed. 

Now  the  single  blade  row  analysis  of  the  previous  section 
demonstrated  that,  for  a blade  row  with  prescribed  characteristics, 
the  flow  distortion  downstream  was  determined  by  the  flow  distor- 
tion entering  from  upstream  along  with  the  knowledge  of  whether  or 
not  the  modified  inlet  flow  stalled  a section  of  the  blade  row  in 
question.  The  relationship  of  the  distortion  in  the  discharge  to  the 
conditions  ahead  of  the  blade  row  is  given  by  expressions  of  the 
type  shown  in  the  third  and  fourth  of  equations  49  and  the  condition 
for  stall  is  determined  by  considerations  similar  to  those  of  equa- 
tions 51  and  52.  There  is  no  difficulty  in  generalizing  this  scheme 
to  one  that  will  apply  to  any  blade  row  of  a multistage  compressor. 

The  generalized  calculation  is  carried  out  in  the  following 
tVi 

steps.  First  suppose  the  n blade  row  is  stalled  over  the  sections 
-ixn/i  6 i . The  coefficients  of  the  downstream  perturbation 

and  the  perturbation  to  inlet  angle  in  the  stalled  region  may  be 
written  in  the  matrix  form 


ft- II  I, 
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The  nine  coefficients  of  this  matrix  are  given  in  the  Appendix. 

Second,  determine  whether  the  value  of  (3„  is  greater  or  less 
than  the  known  value  of  (3*  . J-f  (3*  - ( ^*>  ) then 

the  calculation  above  is  correct  and  the  resulting  information  per- 
mits proceeding  to  the  next  blade  row.  If,  however,  <■ 
the  above  calculation  must  be  repeated  with  the  element 
replaced  by  zero.  Then  this  new  information  permits  proceeding 
to  the  next  blade  row.  The  loss  in  pressure  rise  across  the  stage, 
it  will  be  recalled,  is*  or  zero  depending  upon  whether  the 

blade  row  is  stalled,  and  these  increments  will  accumulate  as  we 
pass  through  the  compressor. 

One  further  point  should  be  mentioned  in  order  to  clarify 
the  application  of  this  procedure;  the  circumferential  coordinate  of 
each  blade  row  has  been  adjusted  so  that  the  distorted  region  appears 
symmetrically  located  with  respect  to  the  local  -axis.  To  recover 
the  physical  flow  field,  the  resulting  profiles  and  angles  should  be 
applied  with  respect  to  an  undisturbed  streamline  passing  through 
the  machine. 


6.  CONCLUDING  REMARKS 


The  example  that  we  have  examined  in  some  detail  has 
demonstrated  that,  within  the  small  perturbation  theory,  the 
exact  solution  consists  of  a linear  combination  of  two  basic  flow 
fields.  In  particular,  the  axial  and  tangential  velocity  profiles 
constitute  a two-parameter  family  where  the  parameters  are 
functions  of  axial  position  through  the  compressor. 

It  is  clear  then,  that  if  we  chose  to  represent  the  solution 
in  terms  of  this  two- paramete r family  for  an  integral  method,  and 
employ  the  actuator  surface  characteristics  of  the  preceding  anal- 
ysis, the  integral  method  will  yield  the  same  solution.  For  example, 
we  represent  the  axial  velocity  in  the  form 
U = L , f u.j)  . M,  l G t*.j  > -•<  ) 

-*  -L,  F C o,  x - $ -i  vx.^G  (<?-,  hv»©(  ) - m J 

where  the  parameters  L,}  r-f,  are  unknown  but  JL,  , are  given. 

The  quantities  'r,/o  and  and  the  corresponding  quantities  down- 

stream of  the  blade  row  have  corresponding  representations.  Inte- 
gration of  the  equations  of  continuity  and  motion  in  the  peripheral 
direction  establishes  the  relationship  between  the  parameters  that 
occur  in  u.,  g-  and -Jo  ; detailed  matching  across  the  actuator  surface 
establishes  the  relationship  between  the  parameters  upstream  and 
downstream.  Note  that  the  matching  across  the  actuator  surface 
corresponds  to  applying  a geometric  constraint  representing  a blade 
shape  or  to  a force  field  representing  a distributed  blade  loading. 

As  a result,  the  differential  equations  that  describe  the  variation 
of  parameters  along  the  compressor  when  the  load  is  continuously 
distributed,  degenerate  to  difference  equations  for  the  compressor 
represented  by  a succession  of  actuator  surfaces.  These  difference 
equations  correspond  to  the  matrix  relations  given  in  equations  48, 
or  more  precisely,  equations  55. 
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APPENDIX 


The  determinant  of  the  matrix  on  the  left-hand  side  of 
equation  48  is 
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The  coefficients  in  the  equations  49  are 
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Figure  2.  Pressure  rise  ^ s a function  of  inlet  angle. 


Figure  4.  Pattern  of  stalled  and  unstalled  regions  on  actuating  line. 
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